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ABSTRACT 

We find the Green’s functions for the accretion disk with the fixed outer radius and time-independent viscos¬ 
ity. With the Green’s functions, a viscous evolution of the disk with any initial conditions can be described. 
Two types of the inner boundary conditions are considered; the zero stress tensor and the zero accretion rate. 
The variable mass inflow at the outer radius can also be included. The well-known exponential decline of the 
accretion rate is a part of the solution with the inner zero stress tensor. The solution with the zero central ac¬ 
cretion rate is applicable to the disks around stars with the magnetosphere’s boundary exceeding the corotation 
radius. Using the solution, the viscous evolution of disks in some binary systems can be studied. We apply the 
solution with zero inner stress tensor to outbursts of short-period X-ray transients during the time around the 
peak. It is found that for the Kramers’ regime of opacity and the initial surface density proportional to the ra¬ 
dius, the rise time to the peak is f^se ~ 0.15 r^^j/vout and the e-folding time of the decay is fexp ~ 0.45 r^y^jvout- 
Comparison to non-stationary a-disks shows that both models with the same value of viscosity at the outer 
radius produce similar behaviour on the viscous time-scale. For six bursts in X-ray novae, which exhibit fast- 
rise-exponential-decay and are fitted by the model, we And a way to restrict the turbulent parameter a. 

Subject headings: accretion, accretion disks - binaries: general - methods: analytical 


1. INTRODUCTION 

The accretion processes are recognized to power many as¬ 
trophysical sources, which are observable due to the release of 
energy when matter is spiralling down into the gravity well. 
The study of viscous accretion disks as underlying mechaisms 
for the extra ction of potential gr a vitational energy began with 
the works oflLvnden-Befll(|1969h : [^akura & SunvaeM(ll973h . 
iNovikov & Thorne! (1 1973L pp.343-450). 

Time-dependent problems for the accretion disks arise 
when describing a wide variety of outburst phenomena, ob¬ 
served in the accreting sources, i.e. binary systems with mass 
transfer and galactic nuclei, and systems with planet forma¬ 
tion. It is usually possible to consider separately the vertical 
and radial structure of an accretion disk due to significantly 
different characteristic time-scales. In such a case the time- 
dependent radial structure of a disk can be described by a 
second-order partial differential equation, which is a conse¬ 
quence of the angular momentum and mass conservation. The 
energy conservation in a geometrically thin optically thick 
disk is provided by the local balance of the viscous heating 
and radiative cooling. 

The basic equation of the viscous evolution, Eq. (ID or dU 
below, describes the viscous angular momentum flux along 
the disk radius. This is a diffusion type equation with a 
variable ‘diffusion’ coefficient. Choice of an analytic ap¬ 
proach to solve the equation depends on the nature of the 
viscosity involved. The preferential way of describing vis¬ 
cosity in astr ophysical disks , although not exclusive, is to use 
a-viscosity (IShakuralll973l: IShakura & Sunvae^ll973h . i.e., 
the proportionality of the viscous stress tensor to the total 
pressure in the disk. One can generally assume that the kine¬ 
matic viscosity v is a product of two power functions: one 
function of radius and another of some local physical param¬ 
eter, the latter usually represented by the surface density. 

The boundary conditions imposed on the accretion disk are 
important. While the outer boundary can be a freely expand¬ 
ing surface due to the spreading of the matter corresponding 


to the outward viscous angular momentum flux, in some sit¬ 
uations certain conditions should be posed at some distance 
from the accreting object. Disks in the binary systems are the 
main focus for such models. 

In a binary system, the tidal torques, acting inside the Roche 
lobe of the accreting object, tru ncate the disk and provide a 
sink for the angular momentum. iPapaloizou & Pringlel (119771) 
show that the effective tidal radius is close to the Roche lobe: 
~ 0.9 of the mean radius of the Roche lobe, well coincid¬ 
ing with the largest non-inter secting periodic orbits in the 
restricted three-body problem dPaczvnskilll977^ . Numerical 
models confirm that most of the tidal torque is applied in a 
narrow region at the edge of the disk, where perturbat ions be¬ 
come nonlinear and strong spiral shocks appear (see [ Pringlel 
Il991tlichikawa & Osakilll994 iHameurv & Lasotall2005L and 

references therein). In the present study, we leave aside the 
outer disk radius variations and other manifestations of the 
tidal interactions in the binary except for the disk truncation. 
We confine ourselves to searching for a solution to a problem 
of the viscous evolution of an accretion disk with the constant 
outer radius. 

We investigate a geometrically thin, axisymmetric Keple- 
rian finite accretion disk with a viscosity in the form v oc r*, 
which is time-independent. We obtain Green’s functions (fT9l l. 
(l22li . (lul l, and (l32l i. which are the kernels in the integrals 
(l20l i. (ITH . and (ITST i to compute the viscous angular momen¬ 
tum fluju F and the accretion rate M. The method of Green’s 
functions allow one to compute E(r, f) and M{r, t) for arbitrary 
initial surface density distributions. Different boundary con¬ 
ditions are considered; no stress or no accretion at the inner 
boundary, which is located at the zero radial coorinate, zero 
or time-dependent mass inflow at the outer boundary. 

It is beleived that the mechanism of bursts in X-ray novae 
involves the thermal instability connected with opacity varia¬ 
tions in the disks, as in dwarf novae. When a burst is ignited, 
and a sufficient portion of the disk is in the hot state, the vis- 

' The couple g of ILvnden-Bell & Pringlel <19740 . 
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cous diffusion starts to govern the evolution, driving the char¬ 
acteristic rise and drop of the accretion rate onto the central 
object. It is a question, which probably allow no universal 
answer, why X-ray novae light curves have different profiles 
and pecularities. If the disk is not very large, as it happens 
in short-period X-ray transients, its whole body is likely to 
be engaged in the viscous evolution. It has been shown that 
an episodic mass input at the outer radius causes the fast-rise- 
expo nential-decay (FR ED) profile due to the viscous evolu¬ 
tion (IWood et al.lboOlh . Using the Green’s function, one can 
explore the situation with an arbitrary initial distribution of 
the matter. Actually, as we show in the present study after 
obtaining the Green function, for realistic initial density dis¬ 
tributions and for the case with the zero inner viscous stress, 
the evolution is generally a FRED. It is interesting to compare 
the results of the model with what happens in the a-disks, the 
prevalent model for the disks’ turbulence in the stellar bina¬ 
ries. 

The second case, when the viscous stress is non-zero at the 
center, applies to the disks with the central source of the an¬ 
gular momentum, e.g., rotating stars with magnetosphers. If 
the magnetosphere’s radius is large enough, the matter has 
a super-Keplerian angular momentum at its boundary. Differ¬ 
ent regimes for such disks have been proposed; propellers and 
‘dead’ disks. Using the Green founction, the disks of constant 
mass that evolve to dead disks can be described. 

In binary systems the mass transfer rate is a key factor. At 
some stages it can be considered negligible comparing to the 
accretion rate in the disk. In other systems, mass transfer vari¬ 
ations may cause the activity in the disk. The case of time- 
dependent mass transfer, which includes the constant transfer 
rate as a particular case, can be also studied analytically if 
V oc r*. Using the Green’s functions, it is possible to express 
2(r, t) and M{r, t) for any outer boundary conditions (for any 
Moutit)) as well as for any initial conditions. 

In Section 2, we present the viscous evolution equation and 
review its known solutions for accretion disks. In Section [3l 
the Green’s functions are found for different inner and outer 
boundary conditions. An analytic way to calculate the disk 
evolution with the variable mass transfer at the outer radius is 
presented; the features of a radiating dead disk are considered. 
In Section @1 a stage of decaying accretion is dealt with and 
a comparison to the cr-viscosity models is made. Application 
to the FRED light curves of short-period X-ray transients is 
presented in Section |5] The last sections are dedicated to the 
discussion and summary. 


by the kinematic viscosity v as follows: 

Wr^^-pvr — . ( 2 ) 

For the Keplerian disks dto/dr - -{3l2)(jjf^lr, and we rewrite 


Wr^{r, t) -3 u)^ I vpdz . 


Zo 

/■ 


With the formula for the surface density 


nr, 


zo 

U) = 2 p(r,z. 


t)dz. 


assuming that v is independent on z, we have 


WrJr, t) - - Wk V £ . 


( 3 ) 


( 4 ) 


( 5 ) 


It is convenient to introduce a new independent variable 
h{r) - v^{r) r - - the specific angular momentum. Then, 

for the Keplerian disk, equation ([T]l can be rewritten as: 

3 (GMf dH^vh) 

dt 4 dh^ ’ h-hK. (6) 


The acknowledged approach is to solve the disk evolu¬ 
tion equation for the variable proportional to the full vis¬ 
cous torque acting between the adjacent rings in the disk, 
cjrip • z 27rr ■ r. This way, the boundary conditions are easy to 
express. From the above expressions, the following relation 
between F -2n Wnpr^ and the accretion rate can be obtained; 


dF 

^n'Lvrr = M{r,t) = — , 
on 


( 7 ) 


where the radial velocity Vr is negative. 

One chooses a procedure to solve (|6]l provided the form of 
V - v(r, 2) is known. This form is determined by the physical 
conditions in the disk. For example, the vertical structure can 
be solved for the a-disk, and the required relation between v 
and S can be found. Let us give a perspective of some studies 
(mainly analytic) of the viscously evolving accretion disks. 
They differ in respect of the viscosity prescription and bound¬ 
ary conditions. The analytic solutions are the benchmarks for 
more sophisticated and involved numerical models. 


2. VISCOUS ACCRETION DISK EQUATION 

The equation of the viscous evol ution in the accre tion disk 
is the diffusion-type equation (e.g.. lKato et al.in~998l) : 


52 _ 1 5 r 1 d 2. 

dt r dr d(tor^)ldr dr ^ 


( 1 ) 


where, among the standard designations of the time, radius, 
angular speed, and the surface density we have the compo¬ 
nent of the viscous stress tensor w^p, integrated over the full 
thickness of the disk 


Wr^r, 0 = 2 



where zo is the half-thickness of the disk at radius r, and the 
stress tensor is related to the differential rotation of gas masses 


2.1. Viscosity v oc r* 

For the kinematic viscosity in the form v oc r*, the equation 
of the viscous e volution becomes a linear partial differen¬ 
tial equation. In ll952lLu^ obtained particular solu tions t o the 
equation proposed by his teacher von Weizsacker (Il948h and 
described the principles to determine the general solution for 
the infinite and finite problems. 

For the accretion disks that can extend infinitely, 
Lynden-Bell & Pringlel ([1973) (hereafter LP74) find, for the 
two types of the inner boundary conditions, the Green’s func¬ 
tions that provide means to describe the whole development 
of the disk. The inner radius of the disk is zero in their exact 
solution. The long-term self-similar evolution of the accretion 
rate goes as a power-law M oc with a parameter I < 1 

(s ee also Append ix A). 

I Pringlel (1199 11) considers the properties of an infinite accre¬ 
tion disk with a central source of the angular momentum using 
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the Gre e n-fun ction technique. A related problem is solved by 
iTanakal (1201 ll) who finds the Green’s function for an infinite 
disk with the boundary conditions posed at the finite inner ra- 
dius. _ 

I King & Ritteil (119981) consider z. finite disk with the viscos¬ 
ity constant over the radius and in time and derive the expo- 
n ential decay of t h e accr etion rate in the disk. 

IZdziarski et al.l (l2009l) study the mass flow rate through a 
disk resulting from a varying mass-supply rate, derive the 
Green’s function for the accretion-rate and its Fourier trans¬ 
form. By means of numerical simulations, they investigate 
a finite disk, residing between non-zero rin a nd rmif, with 
the no n-zero accretion rate at the outer radius. iTanaka et alJ 
(1201 2|) find general Green’s function for an explicit depen¬ 
dence of the inner boundary rin on time and a non-zero mass 
across rin. 

In the present work, we study analytically a finite accretion 
disk with the zero rin and find an analytic solution describing 
its whole evolution, from the rise to decay. For the case of the 
zero stress tensor in the disk center and the mass deposition 
at the ou ter disk edg e , the G reen’s function of the problem is 
found in IWood et all (1200 Ih . 


time-dependent solution with the conserved total angular mo¬ 
mentum. 


3. EVOLUTION OF THE EINITE DISK WITH STEADY 
VISCOSITY 

If the kinematic viscosity v is not a function of E and is a 
function of the radius alone, we arrive at the linear differential 
equation (|6ll. We write the kinematic viscosity as 

V = Vo / . (8) 

Then Eq. (|6]) for the dependent variable F is as follows 


— 

dt 




d^F 

dh^ 


where 

F -2ti — 3 TT/i S vo r* . 

We rewrite (|9]l in the form similar to that in LP74; 

dh^ 4 W / dt ’ 


( 9 ) 

( 10 ) 


( 11 ) 


2.2. Viscosity v oc 

Similarity soluti ons of the first ki nd of the non-linear dif¬ 
ferential equation (lBarenblatllll996h have been found at the 
’’decaying-accretion stage,” when the total angular mo men- 
tum of an infinite disk is constant (lPringlelll974ll98Th . Ac¬ 
cretion rate is found to decrease as oc for the Thomson 
opacity , and as oc for the Kr a mers’ opacity (see also 

lFilipoyiri984t ICannizzo et^ll990l) . iLyubarskij & Shakural 
(119871) obtain self-similar solutions that describe separately 
three stages of the disk eyolution. The first two are the self¬ 
similar solutions of the second kind: the initial moyement of 
the inner edge of a ring of matter toward the center and the rise 
of the accretion rate in the center (M oc for the Kramers’ 
opacity and oc f* for the Thomson opacity). The third stage 
is the decay of the accretion and corres ponds to the solution of 
iPrindel ([1973). The first two stages of iLyubarskii & Shakural 
1198^ may be applicable to the finite disks when the condi¬ 
tions at the distant outer boundary do not affect the behayiour 
o f the disk near the cen ter. 

iLin & Pringle! (!l987l) deriye an analytic self-similar solu¬ 
tion for a disk subject to the grayitational instability, which 
transfers mass and angular momentum in the disk, for which 
M oc and a numerica l solution for the whole eyolution. 
iLin & Bodenheimierl (Il982h found a self-similar solution for 
V oc Y?, if the yiscosity is due to the action of conyectiyely 
driyen turbulent yiscous stresses (M oc 

Accretion disks with a central source of angular momen¬ 
tum and no m ass tran s fer th rough the inner boundary were 
considered by iPringlel (11991 1) who de riyed a self-sirn ilar so¬ 
lution (see also llyanoy et al.l (Il999l) ). iRafikoyl (120131) in the 
extensiye study of the circumbinary disks around supermas- 
siye black holes obtained self-similar solutions for the disks 
with a possible mass transfer across the orbit of the secondary 
black hole. 

For a finite disk, the radial and temporal solution is obtained 
by iLipunoya & Shakur^ (l2000l) yielding M cx f for the 
Kramers’ opacity and M oc for the Thomson opacity. 

In the context of the self-similarity solutions, we would like 
to note that fo r an adyection- dominated infinitely expanding 
accretion flow. iQgilyiel (1 1 9991) obtains by similarity methods a 


where the constants are defined as follows: 

2_ 16/^ 

21 ’ 3vo(GM)F2'' 


( 12 ) 


There is a degeneracy of the index / at = 2. But it is only 
apparent as I enters equation (fTTT i in the denominators. Value 
b - 2 corresponds to the change of the class of functions, 
which satisfy the differential equation, and is left out of the 
following consideration. In the case of b - 2 the yiscous 
time is constant oyer the disk radii fyis P'Ivot^■ We limit 
ourselyes to the finite positiye yalues of I and the Keplerian 
disks. 

In the astrophysical hot accretion disks, the turbul ent yis¬ 
cosity is at work. According to the Prandtl concept (iPrandtll 
Il92-5h . the kinematic turbulent yiscosity is an ayeraged prod¬ 
uct of the two random yalues: Vt = k Ut, where li and Ut are 
the path and speed of the turbulent motion. For the compo¬ 
nent of the yiscous stress tensor we write w^p - -pv^r dcj/dr. 
Following the Prandtl consideration, = -lirdoj/dr, which 

leads to Wnp - pv^. iShakural (119731) proposed the a-disks: 
Wrp - apvl. From the hydrostatic balance one deriyes the 
sound speed ^ or (zo/0 and obtains for a Keplerian disk 
Vt ~ 213) ah{zQl rp. If the relative half-thickness of the a- 
disk is invariable then b - \I2 and I - 1/3. This value of b 
ta kes place for the ‘jS-v iscosity disks’ with v - fih, suggested 
bv iPuschl et alJ (I2000h to act in the protoplanetary and galac¬ 
tic disks. The standard model of the a-disk implies that hjr 
varies with r and Y. 

Let us find the Green’s function of the linear equation (fTTT ) 
with the specific boundary conditions at the boundaries. A 
Green’s function or a diffusion kernel is a reaction of the sys¬ 
tem to the delta impulse, in other words, the Green’s function 
is the solution to the differential equation with the initial con¬ 
dition F{x, 0) = 6{x - xi), where 6{x - .ri) is the Dirac delta 
function. We search a solution by separating the variables, 
F{h, t) - fQioaifi) exp(-s f), where s is a constant and /zout is 
the outer disk radius. Substituting this into (fTTT i. we get 


dV(^) 1 

d^2 4 





(13) 
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where s. This is a Lommel’s transformation of 

the Bessel equation, which can be derived by introducing a 
new independent variable x = A particular solution is 

represented as 

Fii{x, t) - e~^‘ (kx)' [A Ji{kx) + B J^i{kx)] 

(lLustlll95^ . where 7/ and 7_/ are the Bessel functions of a 
non-integer order. If I is an integer number, the Bessel func¬ 
tions of the second kind should be used instead of 7_;. For 
b - 2, the solution is no m ore a Bessel function (see Ap¬ 
pendix I of iKato et al.lll998h . The value b - 2 marks a sit¬ 
uation when the characteristic viscosity time r^/v is constant 
over the radius. We do not consider this case in the present 
study. 

A general solution for an infinite problem is a superposition 
of particular solutions and is expressed by LP74 as an integral 

poo 

over all positive values of k: F(x, t) — Fk{x, t) dk. For a 
problem inside a finite interval, a general solution is a super- 
posit ion of parti cular solutions and is expressed as an infinite 
sum (ILiisllI 19521) : 

CO 

Fix, f) = ^ e-‘'^ ' (ki x)' [Ai Jiiki x) + Bi J^iiki v)], 

(14) 

where constants ki. A,, and B, are to be defined from two 
boundary conditions and one initial condition. The dimen¬ 
sionless coordinate x lies in the interval [0,1], while the phys¬ 
ical coordinate h (specific angular momentum) lies inside 
[0,/Jout]- The outer boundary condition, which sets the ac¬ 
cretion rate at the outer radius of the disk, can be written as 
follows 

dF 

— 4/out(0 at /z — /zout ■ (15) 

ah 

In the simplest case, the homogeneous outer-boundary con¬ 
dition corresponds to the zero mass-inflow rate at the outer 
boundary. As this happens, the viscous angular momentum 
flux F takes some non-zero value, which implies withdrawal 
of the angular momentum at the outer boundary. This corre¬ 
sponds to the situation of the disk in the binary system, when 
the angular momentum is removed from the outer edge of the 
disk by the tidal action of the secondary. 

At .r = 0, one trivial inner boundary condition Fix, t) -0 
approximates the structure for the so-called ‘standard’ disk, 
for which the absence of the viscous stresses at the disk inner 
edge is assumed. Another simple case, dFjdh = 0 at .r = 
0, describes the situation of no accretion through the inner 
disk edge. This case is applicable for the stars with strong 
magnetosphere, when the Alfven radius is greater than the 
corotation radius, impeding the accretion of the matter. We 
consider these two cases separately. 

3.1. The case of zero viscous stress at the inner boundary 
and zero accretion rate at rout 

First, we consider the case of the zero viscous stress at the 
inner boundary. Approximately, this corresponds to the solu¬ 
tion Fix, f) = 0 at X = 0 and only the coefficients at 7;-terms 
in (O remain non-zero. Thus, sum (fT4l l at f = 0 becomes 

C» 

Fix,0)^Yu^kix)‘AiJ,ikix). (16) 

1=1 

We assume the zero accretion rate at rout- Such boundary 
condition approximates the situation of the low mass-transfer 



h/h 


out 


Fig. 1. — Green’s function for the finite accretion disk with the zero viscous 
stress in the center at five moments of time = 0.001, t 2 = 0.01, = 

3/64, ?4 = 0.1, ts = 0.3. Initial spike was at position Xs = = 0.5. 

Constants k = \,1 = 1/3. 


from the secondary in the binary system, for example, when 
the accretion rate in the disk is much higher than the mass 
transfer rate during an outburst. Dropping the time part of the 
solution, we arrive at condition (fTSb required for every term 
in the sum (fThl l. which can be rewritten as follows: 


lJiiki) + kiJ\iki)^0. (17) 


Thus, the general solution for the case of zero viscous stress 
tensor at the center and zero accretion rate at the outer bound¬ 
ary is the sum (fl4l l. where B, = 0 and k, are the positive roots 
of the transcendental equation (fTTl i. To define A, from the 
initial condition, the finite Hankel transforms are used. 

OO _ 

The series 2 k\ Ai Jiiki x) with condition (fTTl l are known 

as the Dini’s series (see IWatsonI 11944 §18.11). Function 
fix) - Fix, 0) x^' can be expressed as the Dini’s expansion 
if the function satisfies the Dirichlet conditions in the inter- 
val and the co e fficients are def ined as kjA,- = 2/y(k;) 7p^(k,) 
(IWatsonlf19441 ISneddonlll951h with the finite Hankel trans¬ 
form 


1 



X fix) Jiiki x) dx , 


0 


where fix) = Fix, 0) x~'. We look for a solution for the 6- 
function as the initial condition: Fix,Q) = 5(x - xi). Using 
the properties of the Dirac function, we obtain 


k\Ai 


= 2x 


i-i 

1 


JiikjXi) 

Jf(ki) 


(18) 


Thus, for the specific boundary conditions, we derive the 
Green’s function, which is the solution to (fTTT i with d-function 
as the initial distribution 


G(x,xi,f) = 2x^x[ ' 


Z 

i 


JiikiXi)Jiikix) 

Jf(ki) 


, (19) 


where k, are the positive roots of Eq. (fTTT i and x = ihjhotffl'^K 
The function is plotted in Fig. [1] for consecutive times. The 
curve for time f 3 corresponds to the maximum accretion rate 
in the center. It is plotted at t-i - given by (1251) below. 
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Fig. 2.— Left: ratio of the time when the accretion rate peaks in the finite disk to that in the infinite disk versus the position of the initial thin ring. Right: ratio 
of the peak accretion rate in the finite disk to that in the infinite disk as a dependence on the initial ring position. 


Given the dimensional initial distribution F{x, 0) in the fi¬ 
nite accretion disk, the distribution of F at any f > 0 is 

1 

F{x, ~ J"G{x, x\,t) dxi. (20) 

0 

The accretion rate can be found from 


x^^ Tout- Applying (l2?t . we obtain 
F{x, 0) = ^ Afdisk 

and the accretion rate evolution is as follows 

M{x, t) - G^(x, Xs,t), 


(24) 


1 

M(x,t) = J T’(xi,0)G^(x,xi,f)dxi//Zout, (21) 


where the Green’s function for the accretion rate 
dG(x, xi, f) 


G^(x,xi,f) = 


dx^‘ 

(xxi)'“^ 






ki 


Jfikd 


( 22 ) 


where we designate 


^vis — 


K" 




21 


T 


r 

'out 

v(rout) 


The time fjnax of the maximum and its accretion rate in the 
disk center Min,peak are of order of those for the infinite disks. 
Using the solution of LP74, we find the time of the accretion 
rate peak for the infinite disk 


4/2 r2 

4(1 +/) “ 3(1 +/) v(rs) 


(25) 


Functions G and G^ are found by iWood et ^ (120011) for the 
case xi = 1 along with an analytical form obtained from the 
asymptotic forms for small and large t. 

One can express the initial F, proportional to the disk mass, 
from the initial distribution of the surface density, using (0 
and (fTSI) : 

16 7r/2 - 

^(2^’0)=^r2 2(r)/. (23) 

where r = h^jGM and h - /zout x^’- 


3.1.1. Spike-like and power-law initial distribution of the viscous 

stresses 


Let the initial mass be located at some radius where the 
Keplerian specific angular momentum is hs = xf ho^t, and the 
mass of the infinitely thin ring is Mjisk, /iout is the maximum 
possible specific angular momentum of the matter. 

The surface density of the d-ring is 


X(r) = 


Mdisk 

2nrs 


S{r- rf) = 


Mdisk 
^TTI r^ 


Xs 6(x - Xj). 


We take into account that 6(r - r^) = 6(x - Xx)dx/dr and r — 


and value (lAlb . The numerical factor in (l25l l is pre¬ 

sented in Table [T]for various, physically motivated values of 
/. In Fig. |2] the ratios of maximum times and accretion rates 
for the finite and infinite disk are shown for several cases of /. 
The ratios approach 1 for Xj « 1 - the disk ‘does not feel’ the 
boundary at the time when the accretion rate peaks. 

If the initial surface density is distributed as a power-law 
with radius L oc we use the following expression in 
Eq. dH); 

Sir ) = («£ + 2) ( r / r 2)°'^'^2 

2n r^ 1 - (ri 

assuming the mass of the disk is enclosed between the radii ri 
and r 2 . 

In Fig.[3]the variation of the accretion rate at the disk center 
is presented for three different initial distributions of F (see 
Fig.lH; two narrow top-hats (a^ = 1/2/ - 5/2) with different 
rj and r 2 and the one corresponding to E(r, 0) oc r. All three 
initial distributions are constructed for the same initial mass 
of the disk Mdisk = 1 ■ The evolution of a thin ring located far 
from the outer disk edge has a characteristic power-law inter¬ 
val, when the disk behaves as if it was infinite. The interval 
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Fig. 3.— The inner accretion rate variation (the black curves with a peak) for three initial distributions of F{h) shown in Fig.|4] The initial mass of the disk is 
the same in all cases. On the left panel, the cumulative accreted mass is shown by the monotonic (green in the electronic version) curves. The dashed curve has a 
power-law interval, which is clearly seen on the right panel. Constants M^isk = T /^out = 1, a: = 1, / = 1/3, ?vis = 3/2. 


lasts about finax[l - (As/Aout)'^']- If 2(r,0) oc r, the mass of the 
disk at f = 0 is enclosed mainly in the outer parts, because it 
grows as the cube of the radius. As expected, the accretion 
rate’s evolution after the peak depends weakly on the pecu¬ 
liarities of the initial distribution of mass in the disk (compare 
the dotted and solid line in Fig. |3]l, as long as its bulk was 
similarly positioned. 


3.1.2. Disk luminosity 

The viscous heating in the disk (its half) is found from F\ 


Gv 


3 F Wk 
Stt 


(26) 


which can be equated to assuming the local energy bal¬ 
ance. 

In the center of the disk, the quasi-stationary distribution 
develops soon and the accretion rate is constant over radii: 
F(t) - M(t) h. The zone of constant M expands over time. 

The viscous heating (l26l l with the quasi-stationary distribu¬ 
tion of F diverges at the zero coordinate. In the real astro- 
physical disks, the inner radius of the disk is a finite value 
Tin + 0. Note that integrating (l26l l over the disk surface from 
tin to Tout also yields the incorrect result for the total power re¬ 
leased in the disk. (As a remark, the integration of (l26l l gives 


o 

&H 



h 


the correct r esult f or the case of the zero accretion rate at the 
center, see § 13.2. 11 1. 

In fact, setting F(rin) = 0 at the finite radius r-m leads 
to a time-independent correction of F near tm- This cor¬ 
rection is found for the inhnite disks by iTanakal (1201 Ih and 
represents the classical factor (1 - V^n/O- Naturally, the 
correction holds for a finite disk. It follows that the bolo- 
metric luminosity can be calculated via the classical formula 
Tboi = G MMl{2ri„) at times t > fvis(nn/^out)^^* with the 
hrst term in the correction of order of (nn/rout) (fvis/0*^^^ 
for earlier times. For the disks in the stellar binary systems, 
usually, ratio nn/rout << 1, and we are safe to calculate Lboi 
by the classical formula with M - M(Q, t) found from (l2lTi . 

3.2. The case of the disk with no central accretion and zero 
accretion rate at rout 

Condition dFjdh -Qd.th-0 necessitates that all A, = 0 in 
(HI, which becomes for f = 0: 


x~' F(x, 0) = Bo + ^ k' Bi J^iiki x) , (27) 

i'=i 


where k, are the positive roots. In the previous section, the 
zero root contributed nothing to (fTST l. Here the zero root pro¬ 
vides a term like (ko xyj^iiko x), which must contribute some¬ 
thing, because the Bessel function /_/ behaves as cx (ko xy‘ in 
the vicinity of zero. 

The method to find B, for i > 1 is analogous to that 
for A/, since the Dini’s series are convergent for the or¬ 
der of the Bessel func tion > -1 (iBenedek & Panzotigll972l : 
iGuadalupe et al.lll993lR We consider only non-integer 0 < 
/ < 1 and obtain 


k‘ Bi = 2 xl 


_, J-i(kiXi) 

A(kd 


(28) 


It is known about the Dini’s series of the negative order 
that an additional first term is to be added to the sum, if the 
equation defining roots k, has particular coefficient^ Just this 


Fig. 4.— Three initial distributions of F(h) in the disk, whose inner accretion 
rate evolution is shown in Fig. [2] Top-hats are located at 0.94 - 1 (solid line) 
and 0.3 — 0.36 of rout (dashed line) and the wide distribution, con'esponding 
to S(r) oc r, goes in (0.1 - 1) rout (dotted line). The disk mass is the same in 
all cases and equals 1. The outer specific angular momentum /lout = T The 
top-hat distributions are constructed using trs = ”1- Parameter / = 1/3. 


^ T he Hankel integral transform was also proved for the Bessel orders 
> -1 ( ^Betancor & Stempakll2001l) . 

^ For the usual form of writing the boundary condition for the D ini’s se- 
ries, H Jy{z) zJ'v{z)=0, the the additional term appears if + v = 0 iWatsonl 

Ell). 
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h/hout ^/\ut 

Fig. 5.— Green’s function for the finite accretion disk with the zero accretion rate at the center at five moments of time ti = 0.001, t 2 = 0.01, ^3 = 0.1 4 = 0.3, 
/g = 0.8 (left) and t[ = 0.001, T 2 = 0.01, = 0.03 ^4 = 0.1, ts = 0.8 (right). Initial spike is at position = (/i//iout)^‘^^^ = 1 (left) and 0.5 (right). Constants k= 

/=l/3,rvis = 3/2. 


case occurs for the properly rewritten outer boundary condi- 
tion (fTSl l: 

(29) 

The additional term in the Dini’s series is expressed as 

1 

Bo = 2(l-l)x~‘ J z^-‘ F(z, 0) z-‘dz, (30) 

0 

where z is a free variable (IWatsonlll944h . Thus, adopting 
F{z) = d(z - xi), from (fl4ll we obtain the Green’s function 
for the equation of the non-stationary accretion with the zero 
inner accretion rate and the finite outer radius; 

G(x, xi,t) = 2(1 - /) x\^^‘ + 

J-l{kiX\)J-.i(kiX) 

(31) 



where ki are the positive roots of Eq. ( |29] | and x - (/j/Ziout)*^^*- 
Formulas (l20l l and (l2Tl i can be applied without changes. 

This Green function is illustrated in Fig. |5] At the late 
stages, the accretion disks forget all information about the ini¬ 
tial distribution of the viscous stresses. For the disks without 
accretion on to the center, the uniform distribution of F devel¬ 
ops, with its magnitude being proportional to the mass stored 
in the disk. 

We take the derivative of (1311 1 and obtain the Green function 
for the accretion rate: 




-w J^i{kiX\)Ji-i{kix) 

out ^ . 


Jtiiki) 


(32) 

It approaches zero when x —> 0 for any time t. This corre¬ 
sponds to our boundary condition at the center: M - 0. This 
disk, without sink of mass and additional mass supply, emits 
radiation even when it reaches the ‘end’ of its evolution, as 
the viscous heating (l26l l does not stop when the flow of mat¬ 
ter ceases. If Fig.|6]its bolometric luminosity is plotted for the 
two initial positions of the ring (the details are given in the 
next subsection). 


3.2.1. The dead disk 


The stable configuration of th e confined ‘dead’ disk (we 
adopt the name from the work bv ISvunvaev & Shakurall977h 
has the following viscous angular momentum flux at each r. 


Fi = F\,^^ - 


4^(l-0/^out‘^'^disk 


= 2 ( 1-0 


Z^out Zffdi.sk 


(33) 

where (l24ll is used and fvis = 8 1 (3 vo)~* as before. 

The total luminosity of the dead disk is a constant value 
Fi (win - Wout)- It can be obtained by the integration of (l26l l 
over the disk surface between rin and rout- Unlikely to the 
case of the zero torque at r = r\n, the late-time value of F 
does not depend on the location at which condition M = 0 is 
set and solely depends on the evolut ion at th e large radii. This 
was shown for the infinite disks by iTanak^ (1201 Ih . Value of 
Hn affects the limunosity only by a change in the size of the 
emitting su rface. There i s still a time-dependent correction 
obtained bv lTanakal(l2011h having an effect at early times t < 

Zvis(On/^out) 


The dead-disk state has to end sometime, as the disk can¬ 
not receive infinitely the angular momentum from the central 
object. If the central star has a magnetosphere, its radius rmag 
can be estimated from the equality of the viscous and mag¬ 
netic torque at the magnetosphere boundary (LP74): 


F(_x 


9 

0) = Kt-^ , 

^mag 


(34) 


where /z* is the magnetic dipole of the statj Kt is a dimen¬ 
sionless factor of order of unity (lLit)unovlll992h . Within 
some factor, this condition is equivalent to the equality of 
the gas and magnetic pressure at the inner disk bound¬ 
ary (ISvunvaev & Shakuralll97'^ . The accretion is inhibited 
if the inner disk radius, rmag, is greater than the corotation ra¬ 
dius, because the drag exerted by the magnetic held is super- 
Keplerian; 


k 'cor 



1.7x m\'1 cm. 


where M*, w* and P* are the mass, angular speed and period 
of the revolution of the central star. In the last expression the 
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Fig. 6 .— Variation of the luminosity of the disk with zero accretion rate 
through the inner radius normalized to value Ljead = <^in- The disk starts 

from the outer radius (the solid line) or from r = 0.5^^ x rout (the dashed line), 
con'esponding to the Green’s functions in Fig. [5] 


mass is normalized by 1.4 Mg. Thus we obtain for a final 
steady configuration 


^ mag,d 
^cor 



^30 vis, 100 


disk,23 


^1/6 
’ OUt ,0 


The ratio can be rather close to unity. The normalizations 
are 10^^ g for the disk mass, 10^° [Gcm^] for the magnetic 
dipole, 100 days for the viscous time, R© for the outer disk ra¬ 
dius. S ituation of > rcor is studied bv iD’Angelo & Spruid 
(I 2 OIOI) . who argue that no considerable expulsion of matter 
from the disk is expected in such regime, as contrasted to the 
propeller scenario. Using the magnetosphere radius as the in¬ 
ner disk radius, we find the steady disk luminosity; 


Ldead^ 1.8X10^^ 


disk,23 


^5/4 3/4 -3/2 „-l^-l/2 -1 

^*^1.4 'out,0 vis,100 ”30 ^ 


for Hn ^ Tout from the both sides, for typical 1 ~ 1 /3 (see 
Table [U. The upper limit on the luminosity of the disk with 
conserved mass is estimated when the inner radius equals the 
corotation radius: 


Tdead,max ~ 4.4X 10^V30^*^'«1.4 ^rg S 

The regime can be sustained for sufficiently strong magnetic 
dipole and fast rotation of the neutron star; 


mP.' ^ 0.63 . 

The flux of the angular momentum F at small radii and the 
inner radius of the disk always manage to adjust themselves 
to the outer variations that proceed on longer viscous times. 
The steady value oc vq could be low due to a degra¬ 
dation of Vq in the outer disk with low luminosity. Then the 
normalized viscous time in the formulae above may be greater 
accordingly to fvis vq. 

Luminosity during the spreading of the ring (Fig. 0 is cal¬ 
culated by integrating expression (l26l l over both surfaces of 
the disk confined between the outer and inner radius; it is plot¬ 
ted in Fig.|6] The inner radius of the disk is varying according 
to Eq. ( l34l i and depends on the star’s magnetic dipole. Its typ¬ 
ical value Tmag.d = Aor ~ 2.4 X 10^^ mZ Rq- 

The torque transfers the angular momentum of the central 
star to the disk and, via the disk, to the orbital motion of the 
binary. The corotation radius increases gradually and, even¬ 
tually , accretion on the center b egins (ISvunvaev & Shakural 
119771: [D’Angelo & Spruit 1201 Ih . The characteristic brake 
time is the angular momentum of the star 27r/*/F* divided 


by the torque F^: 


tbrake 


10^ 


'★,45 


M,'/.Vi{LF*Mdisk,23 


*1.4 out,©-* 


where /*,45 = IflQ'^^ g cm^ and a typical value of I is substi¬ 
tuted. If the mass of the disk grows due to the matter income, 
the braking time of the star shortens, while the magnetosphere 
radius decreases at the same time. 


3.3. The case of an arbitrary accretion rate at rout 

For a variable mass inflow Mout(0 at the outer disk edge 
Tout the solution can be found by a procedure described in the 
Appendix B 


F(x, 0 = J 
0 

/zout rr 

"-JJ 


T'o(-^l) 2 ^out^out(0 


G(x, Xu 0 dxi + 


0 0 


8 //r Mout(T) 
4 U/ hZ 


■ iifout('r) 


X[ G(x, xi, t - T)dxi dr + 


x4/ 

ilfout(0 hout . 


(35) 


where G(x, xi, t) is given by ( fT9l l or OTT l. Substituting G by 
Gyf in the integrals above, and dividing the result by hout, we 
And the accretion rate in the disk M(x, t). We have tried the 
formula obtained in the numerous tests. In particular, it suc¬ 
cessfully reproduces the re sults obtained by numerical simu¬ 
lations of iMineshigd (1 1 994h . The further details, however, are 
out of the scope of the present work. 


4. THE ACCRETION DISK AT THE DECAY STAGE 

We return now to the standard accretion disks with zero 
stress tensor at the center. Eor large f, the first term of the 
sum (1221) dominates, and a single exponential dependence on 
time prevails: 



‘ Jijki xi) 
2im jfikt) 


tk\ 

211 ■ 
^ t fvis 


\ 
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This corresponds to the well-known exponential dependence 
on time of the accretion rate in the disk with the time- 
independent viscosity v. The characteristic time of the ex¬ 
ponential decay is 


t = h^'‘ — 

'exp - “out ,2 


16 fi r^ 

'*4 t 'out 

3 k^ Vout 


(36) 


where we use Vout = 4'orZ- In Table [T] we give locations ki 
of the first zeroes of equation ( fTTl l for several values of 1. The 
table also lists the numerical factors in the expressions for the 
rise and decay time, (|25l) and (O. 

At late times, the profile of F{x, t) is self-similar, because 
G{x,xi,t ^ oo) —> x‘ Ji(ki,x). The distribution E(r) is also 
self-similar. Actually, this profile develops already at the 
very beginning of the exponential decay. This is illustrated 
in Eig. |7]by the fact that the ratio of the inner accretion rate 
to the disk mass, which is the integral of the surface den¬ 
sity, approaches a constant value that depends only on the 
form of the radial distribution. In Table [T] the values of the 
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Fig. 7.— Ratio of the disk mass to the inner accretion rate vs. time normal¬ 
ized by the peak time. The constant value indicates a self-similar distribution 
in the disk has developed. The line styles and parameters are as in Fig. [3] 

parameter aq = /tout/^out can be found. Parameter aq 
de scribes the self-simila r profil e of F{h) and was calculated 
by iLipunova & Shaku^ (120001) for the non-stationary finite 
a-disks. For the stationary disks with non-zero accretion rate, 
Ao = 1. 

What does the presence of this prohle mean for the released 
energy in the disk? Using the energy balance equation for 
the optically thick disks, we obtain that at the outer radius 
= (3/87r) Min i-c. the outer ring luminosity is by 

factor Ao ~ 1.44 less for the viscously evolving disk (/ = 2/5) 
than in the case of the stationary disk. The factor Ao(r) a; 1.16 
for r -0.5x rout- 


4.1. Comparing with the a-disks 

Let us now compare the finite-disk solution for the time- 
independent V and that obtained for the turbulent viscosity in 
the general form v - Vq 2"(f) r'\ For the latter form of viscos¬ 
ity, Eq. (|6]l acquires the following view; 


dF _ d^F 
dt h" dh^ 


(37) 


embracing a dimension constant 


D = 


A -I- 1 
2 


(GMf 


(3 VQ 

\2 {InYiGMf 


l/(a+l) 


and dimensionless parameters m and n: 

A 3 a + 2 — 2b 

m — -, n — -. 

A -b 1 A -(- 1 

For a hnite disk with such type of viscosity the inner accre¬ 
tion rate decays a s a power law (ILipunova & Shakurall2000t 
[Dubus et al.ll200l|) : 

M(f) = M(0)(1+ f/fo)^^^"', (38) 

where M(0) is the accretion rate at the time t - 0, which 
can be attributed to any moment of the decay stage. Eq. (ITtT i 
describes the viscous evolution of the a-disks. 

The characteristic time of the solution to - 
h’'^ l{AmD F'”^^{Q)), where 4 is a numerical constant 
that can be calculated for specihc a and b. Using the 
expression for D and relation Tout = 3 tt hout Tout ^out, we 
obtain 

4 r^ 

to^- -. (39) 

34 A Vout(f = 0) 

Eor the a-disk in the Kramers’ regime o f the opacity, parame- 
ters A = 3/7, b - 15/14, and 4 = 3.137 (ILipunova & Shakural 



Fig. 8. — Relative variation of the accretion rate in two cases: v = vq d’ 
(b = 3/4, the solid line) and v = vq Z" C (Kramers’ opacity, a = 3/1, b = 
15/14, the dashed line). Disk parameters rout and Vout are the same in the 
both models. 


120001) should be used in Eqs. (lJ7ll-(l39ll. In this case, the di¬ 
mensionless factor in (|39]l is close to unity. Value of to de¬ 
pends on the choice of the zero time. 

In the Shakura-Sunyaev disks, the viscous stress is propor¬ 
tional to the total pressure in the dis k, and this prop ortionality 
is parametrized by the a-parameter (IShakurall 19731) : 


Then the kinematic viscosity can be expressed using (|3l as 



Vsound 

Wk 


(40) 


or 


Vt 


2 

-aa)y,r 


2 



(41) 


To approximate the viscous evolution of the a-disk by 
the disk with steady viscosity, one has to estimate the ap¬ 
propriate parameter b for Eq. (I8]l and (|9]l. This can be ef¬ 
fectively done using the relation Vt oc r'^^(zo/r)^ oc r*. 
The solutio n for the stationary press ure-dominated station¬ 
ary regions (IShakura & SunvaevI 19731) gives the relative half¬ 
thickness zo/r oc consequently, b ^ 3/4 or / 2/5 

in (O, for the free-free opacity regime. Table [T] lists the 
parameters approximating in a similar way the standard a- 
dis k with the Thomson opacit y and the ‘dead’ disks mod¬ 
els (ISvunvaev & Shakuralll977l) . 

In Eig. [8] the two laws are plotted; exponential oc 
exp(-f/fexp) for / = 2/5 and power-law (l38l l for m = 3/10. 
Viscosity at the outer radius and the outer radius itself are the 
same. If applied to the real bursts during the time spans of or¬ 
der of ~ fexp, the mod els cannot be easily discriminated (this 
was also pointed at bv iDubus et al.ll200ll) . This fact provides 
a possibility to estimate a using the e-folding time of the light 
curve. In the next Section, we show how it could be done us¬ 
ing the observations of the disks in X-ray novae, whose light 
curves in many cases show exponential decays. 


5. THE BURST LIGHT CURVES OE X-RAY NOVAE 

In the X-ray novae, binary systems with a compact star 
and, typically, a less-than-solar-mass optical star, transient 
outbursts are th ought to orig inate due to a dwarf nova-type 
instability ('e.g.. lLasotall2001i) . The temperature rises in a por¬ 
tion of the ‘cold’ neutral disk, which becomes ionized, and the 
viscosity v rises. Moreover, the turbulent parameter a changes 
from acoid to higher ahot- The heating front moves outward 
from some ignition radius with the speed Ufront ~ (Thot Vsound- 
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TABLE 1 

Parameters of the Green’s functions for two types of inner boundary condition. 


b 

(1) 

1 

(2) 

ki 

(3) 

Cx(''s/''s) ‘ 

(4) 

4xp(^out/^out) 

(5) 

flo 

(6) 

Comment 

(7) 

0 

1/4 

1.0585 

1/15 

0.298 

1.267 

constant v 

1/2 

1/3 

1.2430 

1/9 

0.383 

1.363 

a-disks with constant hjr, /3-disks 

3/5 

5/14 

1.2927 

0.125 

0.407 

1.392 

fT-disks, Thomson opacity (^io = 1.376^) 

3/4 

2/5 

1.3793 

0.152 

0.449 

1.444 

fT-disks, Kramers’ opacity {ao = 1.430^) 

I 

1/2 

1.5708 

2/9 

0.540 

1.571 

F(h) oc sin((;r/2)/i//iout) 

2 

oo 

— 

— 

— 

— 

/vise does not depend on r 

2/5 

5/16 

3.4045 

0.189 


0 

dead cr-disks, Thomson opacity 

3/5 

5/14 

3.3425 

0.265 

— 

0 

dead cr-disks, Kramers’ opacity 


Note. — Columns are as f ollow s: (1) The index of the radial dependence of v; (2) I defined by <12L (3) The 
first zero of equations (Tt) and i29\ . in the upper and lower part, respectively; (4) The dimensionless factor in j25t 
(the upper part) and the one in a similar formula derived by LP74 for the time of m aximum F in the disk without 
central accretion (the lower part of the Table); (5) The dimensionless factor in <361 ; (6) Parameter describing the 
self-similar radial profile, ag = Min/lout (7) Note on the applicability. The opacity law is indicated for the 
ff-disks. The applicability of parameters is approximate for the tr-disk. Refs: iLipunova & Shakural <20001) . 


It stops at a radius, where the surface density of the disk 
is too lo w to sustain t he stable (on the thermal time-scalel 
hot sta te (iMeverll 19841 iMenou et al.ll200Cil: [Dubus et al.ll2001t 

lLasotall2001l) . 

In the dwarf novae, the hot state accretion is quenched 
swiftly (comparing to the viscous time) by the matter return¬ 
ing to the cold state. It was suggested that in the X-ray no¬ 
vae the illumination from the center has a stabilizing effect 
on the outer disk, keeping the temperatures above the hy¬ 
drogen ionization temp e rature d Meyer & Me yer-Hofmeisteij 
11984 IChen et aT]ll993l) . iKing & Ritteif(if998h propose that 
in order for the disk to undergo the viscous evolution an in¬ 
tensive irradiation must be included. 

A canonic al light curve of an X-ray nova is described 
as a FRED (IChen et al.l 119971) . Theoretically, if the disk 
is truncate d from outside, a steeply decaying light curve is 
produced (iKing & Ritted 119981: iLipunova & Shakural 120001: 
IWood et al 120011) . The main property is that the disk does not 
expand infinitely, which is a consequence of its being in a bi¬ 
nary system, and i ts radius does not cha nge on a viscous time 
scale. Scenario bv iKing & Ritterl (119981) produces exponential 
decays due to the steadiness of the parameter v in time. More 
realisti c temporal dependence of v leads to a steep po wer-law 
decay (ILipunova & Shakurall2000l: iDubus et al.ll2001l) . which 
is close to an exponential one during few viscous times (see 
the previous Section). 

The possibility to model the X-ray novae bursts, which 
have the FRED profile, by the viscously evolving disk with 
the steady visco sity and afte r a des crete mass-transfer event, 
w as shown by [Wood et al.1 ( 120011) (see also the modeling 
in ISturner & Shradeiil2005l) . If a burst is due to an instability 
affecting a slab of matter extended along the radii, the Green’s 
function obtained in the present work should be used. Effec¬ 
tively, after a heating front has passed, the dwarf-nova type 
instability leads to some distribution of the hot-ionized matter 
with high viscosity. In the most trivial case, when the bulk 
mass is concentrated near rnnt, the evolution of the disk ap¬ 
proaches the one found bv lWood et al.l (120011) . 

Let us compare the typical time of the heating front prop¬ 
agation r/Vfront and the time f^ax it takes the inner accretion 
rate to peak, which is of order of (see (l25l l and Fig. |2]i. 


Relating a and v by (l40l l. we obtain: 


tfi-ont _ ^ + 1 r^'^Lund lZi_Gi 

fmax 2P y GM ' y nm„ ' 

where the gas temperature is normalized by 10"^ K, central 
star mass by a solar mass, and radius by lO" cm, p is the 
molecular weight of the hot matter. Radius is the charac¬ 
teristic radius where the bulk mass resides at the beginning of 
an outburst. This ratio indicates that we are safe to model an 
X-ray nova outburst at time ~ f^ax using the Green’s function 
obtained. 

If one assumes the commonly suggested surface density 
di stribution in the quiescent state in the X-ray novae, S oc 
r (lLasotall2001l) . then, as it is shown in Fig. 12 after the peak 
of the burst, the only important parameters are the disk mass 
and the outer radius. For such density distribution, the bulk 
mass of the disk resides close to the outer radius. As we show 
in Sect. @1 already at the maximum of an outburst the disk 
acquires the self-similar radial distribution, which is not de¬ 
pendent on the initial distribution of matter. The necessary 
condition for this is that the bulk mass of the viscously evolv¬ 
ing disk is located near the outer radius. This is consistent 
with the absence of a power-law interval before the exponen¬ 
tial decline on the the light curves of the X-ray novae that 
show the FRED behaviour. 

Eqs. (|25] | and (l36] l yield the ratio fexp/fmax taking into ac¬ 
count that fmax is very close to for ~ rout- The value of 
fmax/Cax i^ depicted in Fig. |2 where fmax is the peak time 
for the finite disk. For / = 2/5 (see Table [T]i we obtain 
fexp/fmax ~ 3. This provides a rough test whether the par¬ 
ticular light curve is produced by the viscously evolving disk 
made of hot ionized matter and with approximately constant 
outer radius. 

The viscous-disk relation frise ~ 0.15 for the hot- 

disks with the Kramers’ opacity in the outer parts leads to a 
conclusion that to obtain a secondary maximum on the X-ray 
light curve of an X-ray nova, the disk has to provide an extra 
mass input at the radius that is determined by the time of the 
input. The secondary maximuma are usually observed around 
(1-2) fexp after the peak. If the extra mass input happens at the 
time of the burst maximum then it should happen at least at a 
radius about 2.4 rout (from fexp/frise ~ 3 and v oc r^^'*), without 
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Fig. 9.— The peak-normalized light curves of GRO J0422-(-32 in 1992, A 0620-00 in 1975, GS 1124-68 in 1991, GS 20004-25 in 1998 (as collected in IChen et al.l 
119971) and 4U 1543-47 in 2002 and XTE J1753.5-0127 in 2005 (the quick-look results provided by the ASM/RXTE team). The energy bands of the X-ray data 
are indicated in the plots. The model curves are the peak-normalized accretion rate, calculated by 1211 using t^xp indicated for each burst, and 1 = 2/5. The initial 
surface density distribution E oc r, and the inner radius of the initial hot zone is 0.01 X rout- Note the different axis scale for A 0620-00, for which two different 
models are plotted, with the initial inner radius at 0.001 X r^ut (solid line) and 0.3 X rom (dotted line). 




Fig. 10.— Top panel: estimates for the turbulent parameter a times the 
square relative half-thickness of the disk at rout for the six bursts of X-ray 
novae. Lower panel: corresponding estimates of the viscosity parameter v at 
r =Ro. 

mass input from the intermediate radii. Another possibility is 
that the secondary maximum is triggered later after the peak 
and then it can happ en at a radius ~ rout (as in the model of 
lErtan & Alt)aiil2003) . 

If an X-ray nova has the viscously evolving accretion disk 
with constant rout during the interval of the ‘exponential de¬ 
cay’ (usually, first ~ 50 days after the peak), we can infer the 


value of the turbulent parameter a, using the closeness of the 
model with steady visosity and the a-disk. As the irradiation, 
which provides the steadiness of the hot zone outer radius, 
does not change the vertical structure of the disk (iDubus et alJ 
mn), we can adopt / = 2/5 applicable for the hot ionized 
standard a-disks. Combining Eqs. (l36]) and (HTI) . we get 


O' 0.15 


^out ^ I M \ ' 

2Rej \ 0.05 ) \10Mo/ IsO'i/ 


, (42) 


which relates the a-parameter to the observed fexp, the 
e-folding time of the bolometric luminosity or the inner ac¬ 
cretion rate. The last formula resembles quite a few others 
found in the literature. Similarly, the relation of a to the half¬ 
thickness of the disk cannot be eluded. However, an advatage 
is that (l42l l contains the e-folding time which can be inferred 
quite accurately from observations. 

We have sampled six ERED light curves in different X-ray 
novae (Eig. |9]) and fitted them with the steady-v viscously- 
evolving disk model. Eor particular values of fexp, the model 
curves agree well with the peak-normalized count-rates be¬ 
fore the second maximum, which is observed in most of the 
cases. Eor A 062 0-00, we use the value of fexp, estimated by 
IChen et ^ (119971) from the 3-8 keV light curve during the de¬ 
cline. Variation of the initial inner radius of the hot zone leads 
to slight changes of the model curve for the particular initial 
distribution of 2(r), as can be seen in the panel for A 0620-00. 

In Eig. [To] we plot the estimates for cr-parameter times the 
normalized square of the relative half-thikness at the outer ra¬ 
dius of the disk. Theses estimates are obtained using (l42ll 
with the values fexp indicated in Eig. |9| In the lower panel, 
the viscosity v is plotted, found accordingly to (l36l) . As v is 
not constant over radius, we choose to present its value at the 
same radius R© for all cases. The considered X-ray novae 
have known orbital periods. The masses of t he primary and 
secondary components, mass ratios are from [Charles & Coa 
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( |200^ . For bla ck hole candidate XTE J1753.5-0127 with 
Poib = 3.24 hr (IZurita et al.ll2008h we set rather arbitrarily 
M - 3 - 15 Mq and Msec = 0.2 - 1.2 Mp. The Roche lobes 
/?RL are calculated using the formula bv lEgglet^dl983h . and 
'"out/^RL = 0.88 + 0.02. The vertical barrs in Eig. [TOl reflect 
the errors related to the uncertainties of the binary parameters 
a nd do not take into accou nt uncertainty of fexp. 

ISuleimanov et al.l(l2007h give the half-thickness of the outer 
parts for the stationary a-disks when the hydrogen is com¬ 
pletely ionized; 

zo/r ~ 0.05 (M/Mo)-3/« m3/20 ^i/ 8 ^-i/io ^ 


with the parameters’ normalizations equal to lO'^gs ' and 
10'® cm. It can be rewritten for the factor in Eig.[T0]as follows 


( 


£^y 

0.05/ 


^ 1 X 


M 

lOMoj 


M 

O.lMeJ 



a, 


•1/5 

0.1 ■ 


where a is normalized to 0.1 and Mcr is the conventional rate 
of accretion producing the critical luminosity; 0.1 - 

1.3 X 10^®(M/Mo) ergs^'. Eor the black hole masses in the 
range 3-15 M©, accretion rates in (0.01 - 0.5) Ma, a in 
0.1-0.5, and rout in 1-3 R©, the above factor changes between 
extremes ~ 0.3 and ~ 3.8. Eo r 4U 1543-47 with parameters 
MjM rr ^ 0.25 (IWu et al.ll2010l) . M = 9.4 M© (ICharles & Coel 
120061) . rout = 4.5 R© = 0.9 X Rrl and a - 0.2 one obtains the 
factor (zo/r/0.05)^ » 1.8. It follows from Eig.[T0]that a can 
be less than 1 in this case, but this is still not a self-consistent 
value. Apparently, the disk in this longer-period X-ray nova 
was not entierly involved in the viscous evolution that pro¬ 
duced the ERED light curve in 2002. Levels of a of the other 
five bursts (Eig. [TO] the upper panel) agree with the estimates 
for the fully ionized, time-dependent accretion dis ks obtained 
by various approaches in X-ray and dwarf novae (iKing et al.l 


120071 ISuleimanov et al.ll2008l; IKotko & Lasot^l2012h . in the 
decretion disks of Be stars (ICarciofi et al.ll2012). 

The values of fgxp, indicated in Eig. |9] are to be considered 
with caution. Note that the observed light curves are limited to 
an energy band. Spectral modeling might be required in many 
cases to obtain accurate fgxp. If an X-ray nova during the burst 
is in the high/soft state, when the disk emission dominates the 


spectrum, the spectral modeling can provide Tin(t) oc My^(t) 
in a straightforward way. We expect fgxp not to alter by a factor 
larger than 2 for the X-ray flux dominated by the disk (prelim¬ 
inary resluts for 2-10 keV). 

On the other hand, some X-ray novae bursts proceed en¬ 
tirely i n the low/hard state, su ch as the burst of XTE J1753.5.- 
0127 (iMiller & Rvkofll[2007^ . Still, the form of the corre¬ 
sponding X-ray light curves is informatve. Even if the opti¬ 
cally thick disk does not reach the innermost orbits, it is very 
probable that the radiated energy is proportional to the accre¬ 
tion rate, because the released potential gravitational energy 
is the ultimate source to power the source. Second, the inner¬ 
most disk peculiarities do not influence the behavior of M(f) 
as the latter is defined by the process of the viscous matter 
redistribution in the outer disk. 


6. DISCUSSION 

Equation (|6|l is written without the tidal-stress term, fol¬ 
lowing the proposition that it should be negligible everywhere 
in th e disk, except in a thin ri ng near the tid al truncation ra - 
dius (llchikawa & Osakilll994ft . According to lPringl^ (119911) . 
effects near the tidal-torque radius can be approximated with 


an effective boundary condition M - 0. Numerical models 
confirm that most of the tidal torque is applied in a narrow 
region at the edge of the disk, where perturb ations become 
non-linear and strong spiral shocks appear (iPringlel 119911; 
llchikawa & Osaki[ l994tl Hameurv & Lasotall^05L and refer¬ 
ences therein). As Smakl(l2()0'^ notes, in the case of dwarf no¬ 
vae, during outbursts, the observed values of the outer radius 
appear to be approximately consisten t with the theoretical pre¬ 
dictions. iHameurv & Lasotal (l2005h argue that the action of 
the tidal torques are important also well inside the tidal ra¬ 
dius. They notice at the same time that observations do not al¬ 
low discriminating between the rapidly growing and smoother 
tidal stresses, as the model light curves are not strongly af¬ 
fected. 

The tidal radius for all mass ratios can be approximated as 
0.88 + 0.02 of the average Roche lobe radius of a given com¬ 
ponent in circular binaries (iPapaloizou & Pringlelll977]) . The 
disk can extend be yond this radius due to the high inclination 
of the binary orbit (IOkazakill2007^ . Effects of the high eccen¬ 
tricity of the binary orbit on the disk truncation radius have 
been studied by various methods (e.g..lArtvmowicz & Lubowl 

IT^IPichardo et alJl2()()llOkazakill2007l); 

The gravitational influence of the secondary on the disk in a 
binary system leads to a number of interactions, non-resonant 
(or tidal) and resonant, any of those tr ansferring angular mo¬ 
mentu m from the disk to the binary (lArtvmowicz & Lubowl 
Il994ft . The tidal distortions generally lead to the truncation 
of the disk, whereas the resonances of different strength im¬ 
pact or do not the disk str ucture, depending on thei r growth 
rate a nd other conditions (I Whitehurst & Kirill 19911; iQgilviel 
|2002). An example of a resonance action is the deb ated ex¬ 
planation of the superhamps in SU UMa stars (see, e.g. [Lubowl 
119911; iKornet & Rozvczkdl2000l; IHameurv & Lasotall2005h . 

In our mathematical set-up, we have presumed a fixed outer 
radius of the disk without going into details of its actual value. 
Thus the model is appplicable to the coplanar disks in low- 
eccentricity binaries. The outer radius of the disk can also 
change due to intrinsic processes. Reproducing the whole 
evolut i onary cycles of the dwar f novae and X-ray novae. 
Is maid (l200(]h . iDubus et al.l (120011) should have taken into ac¬ 
count the variations of the outer disk radius, but these changes 
are significant during the time spans much greater than the 
viscous times. 


7. SUMMARY 

Under some astrophysical conditions, viscously evolving 
disks formed in binary systems are effectively truncated from 
outside. In the present work, we use the method proposed by 
lLusti(ll952h to find the Green’s functions for the linear viscous 
evolution equation for the disk with the finite outer radius. 
Green’s functions are obtained for the viscous angular mo¬ 
mentum flux F and the accretion rate M. Two inner boundary 
conditions at the zero coordinate are considered. They cor¬ 
respond to the accreting disks and to the disks with no mass 
sink at the center. The Green’s functions allow one to compute 
E(r, t) and M{r, t) for arbitrary initial surface density profiles 
that develop into self-similar distributions on sufficiently long 
time scale. The analytic formula to calculate the disk evolu¬ 
tion with the variable outer mass inflow is derived. A solution, 
which can be found with the use of the Green’s function, has 
the properties defined by the main equation of the viscous evo¬ 
lution and thus provides the basic description of the transient 
phenomena in the viscous disk. 

The Green’s functions are found in the form of quickly 













































































13 


converging series and can be easily reproduced with standard 
computer methods. The C-code written for the case of the ac¬ 
creting disc with the use of the GNU Scientific LibrarjQ can 
be downloadecQ. 

We present the relations between the rising time, e-folding 
time, and disk viscosity. For six bursts in the X-ray novae, 
which are of FRED type, we show that the model describes 
well the peak-normalized light curves before the second max¬ 
imum. This favors the mainly viscous nature of their evolu¬ 
tion during this period and enables us to obtain an estimate of 
viscosity v, which depends on the outer disk radius. 

The models with time-independent viscosity are shown to 
approximate well the evolution of a-disks during time inter¬ 
vals comparable to the viscous time. Consequently, estimates 
of the a parameter may be derived for the disks in the bi¬ 
naries with known period and masses. These estimates rely 
on the relative half-thickness of the disk at the outer radius. 
The estimates of a for the five short-period X-ray novae out¬ 
bursts, GROJ0422H-32in 1992, A 0620-00 in 1975, GS 1124- 


68 in 1991, GS2000 -h 25 in 1998, and XTE J1753.5-0127 
in 2005, are in lin e with the values estimated so far for the 
hot viscous disks (iKing et^l2007t ISuleimanov et akl 120081 : 
iKotko & Lasotall2012^. 

Another Green’s fucntion is found for the disk that has zero 
accretion rate at the inner radius and acquires angular mo¬ 
mentum from the central star. In the steady state, the mass of 
the disk cannot be very high and the disk has low luminosity. 
It radiates most of the rotational power transferred from the 
central star. This disk works as a gear transmitting the angu¬ 
lar momentum of the central star to the orbital motion. The 
‘dead’ stage is ended by an abrupt fall or trickle of the mat¬ 
ter on to the star after the centrifugal barrier has moved close 
enough or beyond the magnetosphere’s radius. 
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APPENDIX 

INEINITE DISKS: GREEN’S FUNCTION AND BASIC RELATIONS 

Green’s function for the infinite disk, obtained by LP74, can be written in the dimensionless form with our designations (c.f. 
formula (fT9l l) as follows: 

f f / V V \ 

/r, xxj J x,+x ^ 


G°°(x,xi,t) = 


2f 


■ exp 


4f 


where Ii is the modified Bessel function of the first kind, x = is some characteristic value of the specific angular 

momentum, which we use instead of /lout- 

To find the physical distribution of F{x, f) in the course of the evolution of the initial narrow ring of matter, initially located at 
the coordinate Xs, one takes the integral of F).!!, 0) (l24ll with the kernel on the interval xi 6 [0, oo]: 


77 ... _ Mdi^^hcl(xx,)‘ ___ 

^^ 6X.p I 


At 




where is the initial mass of the ring. 

The accretion rate at the center can be found from Min = dFldh\h-,o (see (ISTT l). One finds 


, 1 - 2 / 


M.„(f) = 


dF\ 


MdiskU e 


.1 „-Tjt 


or 


21/ic lx-»o T(l) 
i+z 


+ 1 +/ 


A^m(0-M“peak(Y) 


e 


-Teh 


where the typical times Tpi and Tg are introduced. 


.2 Gh 


h 


1 +l 


Tg = 


■Tpl. 


4 e 

We find that the accretion rate peaks at the time (l25l l - Tgl{\ 1) with the value 

iUdisk (1 + 0^ 


m; 


in, peak 


Cx ei+'F(0 


(Al) 


SOLUTION EOR THE CASE WITH VARIABLE ACCRETION RATE AT THE OUTER EDGE OE THE DISK 

For infinite disks, see lMetzger et alJ (120 1 2h and lShen & Matzri^ (l20I4h . who have obtained the Green’s function solution to 
the viscous evolution of a disk with mass sources/sinks, which are distributed in a fashion along the disk radius. 


http://www.gnu.org/software/gsl/ 


’ http://xray.sai.msu.ru/~galja/lindisk.tgz 
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We have the following problem with the inhomogeneous boundary and initial conditions: 

d^F 1 


dh^ 4 \ l I dt 
dFldh\h=h„„, = MoutiO 
F\,=o^Fo(h). 


(Bl) 


Let us substitute function F(h) by a function F{h) using the relation 


F(h) = F(h) + 




2h„ 


■ 4^out(0 ■ 


(B2) 


This substitution is used, for ex ample, to solve the problem of finding the temperature distribution in a cylinder, on a surface of 
which a thermal flux is defined (iBogolvubov & Kravtsovl[r998l) . For the new F(h) the problem is the following: 


dF 


= 4 




dt 

dF/d!l\h:=:h„, = 0 , 

F\,^o = Fo(h) = Fo(h)- 


(B3) 




MUt = 0), 


2 /Zout 


where 

(bth ft = 4 I - 1 — 

, , ^out 2 /Zout tit 

The corresponding Sturm-Liouville problem with a free variable ^ = /z//zout is investigated in Sect. [2 it consists of Eq. (fOT l 


5^2 




(B4) 


where s,- or kj = plays a role of an eigenvalue, and the boundary condition dfi(^)ld^ = 0 at ^ = fout- Changing to the free 

variable x = ^*^2' = (/z//zout)*^^^ we can write for a non-integer / 

fi{x) = ikix)‘ [A; Ji(ki x) + Bi J-i(ki x)], 

where A, and B, depend on the inner-boundary condition. According to the Steklov theorem, the solution to (IB3l l can be obtained 
in the form: 


F(h, 0 = ^ Ui(t)fi(x ), 


(B5) 


whe re co effici ents Ui{t) are to be found from the inhomogeneous first-order differential equation over t. Let us substitute (IB51 l 
into (IB4b and (IB31 l. multiply by fn(x) and integrate over x from 0 to 1, using the orthogonality property of the eigenfunctions. We 
arrive at 


dun(t) 


dt 
lf=0 “ 


H- Si Unit) = d)„(f) , 


with designations 


= J Foix)fnix)dxl\\fn\\^- <b„(f) = J <i>(h(x),t)f„(x)dx/llf„f , 
0 0 

where ||/„|p is the norm of eigenfunction /„. Solution to (IB61 > is as follows 

t 

Uiit) = 0/ ^ + J' Oi(T) dr. 

0 

Let us write down the solution to (IB3l l for the case Mout = 0, that is, <I)(/z, t) = 0: 

CO 

F(x,t) ^ fix). 


(B6) 


(B7) 


1=1 
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The last expression is equivalent to (fTO . In Sect. |3] we find the coefficients 0, = A, and B, for the Dirac delta function as the 
initial condition and express F{x, t) using the Green’s function. In the analogy, the solution to the problem with Q>{h, t) Q has 
the view 

1 / 1 


/ 


F{x,t) — I Fq(xi) G(x, Xi, t) dxi t) G(x, Xi,t - t) dxi dr 


//■ 

0 0 


(B8) 


where G(x, Xi, f) is given by ( fT9b or ( ISTT i. 

Substituting (IB8b into (IB2b and using functions Fq and d) from (IB3b . we finally arrive at dTSl l. The coordinate integral in its 
second term can be taken analytically; it involves the Lommel functions and can be tabulated beforehand for the particular values 
of I and ki for the sake of the computational efficiency. The second term of dTSl l for B, = 0 can be expressed as: 


where 


/ZnMt X 


I 

^ Ji{ki x) J" 


exp 


t-T k. 




V ^vis 2/ y 


41 Lx —^- L2 Mout(T) 


dr, 


Lx = 


= ' Ji{kiXi)Axi jjfiki)-, ^2 = J" Ji(ki xi)dxi jjfikj) . 


(B9) 
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